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Effect of non-uniform exchange field in ferromagnetic graphene 
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We have presented here the consequences of the non-uniform exchange field on the spin transport 
issues in spin chiral configuration of ferromagnetic graphene. Taking resort to the spin orbit coupling 
(SOC) term and non-uniform exchange coupling term we are successful to express the expression 
of Hall conductivity in terms of the exchange field and SOC parameters through the Kubo formula 
approach. However, for a specific configuration of the exchange parameter we have evaluated the 
Berry curvature of the system. We also have paid attention to the study of SU(2) gauge theory of 
ferromagnetic graphene. The generation of anti damping spin orbit torque in spin chiral magnetic 
graphene is also briefly discussed. 


I. INTRODUCTION 

Graphene [1] is a two dimensional (2D) layer of carbon atoms with unique flat geometry and remarkable intrinsic 
transport properties. Various unusual properties of it make graphene an important material for the study of spin 
physics. The long spin flip length of graphene has made it a potential candidate for spintronic [2, 3] applications 
where the spin orbit coupling (SOC) plays a very important role. 

In recent times, magnetized graphene has attracted a lot of attention from application standpoint. There are variety 
of ways for the experimental realization of magnetized graphene or more precisely graphene with spin imbalance. 
There may exist some intrinsic ferromagnetic correlations in graphene. Use of an insulating ferromagnetic substrate 
or adding a magnetic material or magnetic dopants or defect on top of the graphene sheet may be other options to 
achieve ferro-magnetism in graphene. In particular, by depositing ferromagnetic insulating layer (FI) on graphene, 
magnetization is induced through the exchange proximity interaction (EPI). Induction of large exchange splitting has 
been demonstrated by depositing ferromagnetic insulator EuO on graphene [4]. Thus induced magnetic effects in 
graphene basically rely on the principle of exchange interaction. 

Here we have considered the spin chiral nature of ferromagnetic graphene [5] , where the sub bands of charge carriers 
with opposite spin direction are of electron or hole like. So there exists a regime in which majority and minority spin 
carriers belong to different bands i.e conduction and valence bands. In this spin chiral configuration, the authors in [5] 
demonstrate a different type of Hall effect from the spin Hall mechanism in the extended Drude model approach. The 
Hall effect that is generated by the spin Hall mechanism in graphene via the spin-orbit dynamics of the carriers can 
be well explored through non uniform exchange interaction [6-8]. As the gauge theory can be applied in describing 
transport in graphene [9, 10], the role of momentum space and real space Berry curvature for a system with finite 
spin-orbit coupling is also worth studying. 

On the other hand, in magnetic memory technologies there is a vital role of spin torque in magnetized switching via 
spin Hall effect. Spin transfer torque (STT), an important aspect of spintronics has been studied from the perspective 
of local magnetization dynamics. Recent theoretical and experimental studies suggest that the SOC can induce STT 
unlike the conventional spin injection based STT. Simultaneous presence of SOC and local magnetization will induce 
a SU(2) gauge potential. Very recently, in ferromagnetic semiconductor an antidamping spin orbit torque (SOT) has 
been observed [11] which originates from the Berry curvature and provides an analogous understanding of intrinsic spin 
Hall effect via Berry curvature. Extension of this formalism for ferromagnetic graphene may provide some important 
consequences in the study of spintronic applications. 

In this paper, we would like to study the spin orbital dynamics of the carriers in spin chiral ferromagnetic graphene 
[5, 12] within a framework of semiclassical theory. We find that the non-uniformity of the exchange coupling can be 
a potential candidate to explain various spin related aspects in spin chiral ferromagnetic graphene. 

In case of spin Hall effect (SHE), the SOC provides the gateway to study different spin related issues. The well 
known Kubo transport formula has been used to study the Hall conductivity of spin orbit coupled system. We have 
explicitly used this approach to calculate the conductivity of a ferromagnetic graphene system, where the exchange 
coupling is momentum dependent. In our system the total conductivity arises from the SOC due to external electric 
field and non-uniform exchange field. It has been shown that the Hall conductivity can produce pure spin current in 
our magnetic graphene system. 
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Furthermore, the k dependence of the exchange coupling may induce a different kind of SOC. Due to the presence 
of these two kind of SOC, we have a SU(2) gauge field in the system, from which we can calculate the Berry curvature, 
which is different for the up and down electrons. 

The paper is organized as follows: In sec II, the Hamiltonian for the ferromagnetic graphene is formulated. The 
derivation of Hall conductivity is shown in the Kubo formula approach in section III. In sec IV We have studied the 
spin gauge and curvature for a particular configuration of momentum dependence of the exchange field. Section V 
deals with the derivation of Berry curvature for coordinate dependent exchange coupling. The section VI includes the 
derivation of a anti-damping spin orbit torque, when the exchange field is considered as time dependent. We conclude 
in section VII with a discussion about the importance of the non-uniform exchange coupling. 


II. THE HAMILTONIAN 

We here consider the spin chiral configuration of the ferromagnetic graphene. In this spin chiral configuration 
the sub bands of charge carriers with opposite spin direction are of electron or hole like. When we consider such a 
configuration of ferromagnetic graphene with a substrate induced SOC, the corresponding Hamiltonian can be written 
as [5] 


H = vpd.k + Ep + a.h + V(r) + X g[<x x V r V(r)].fc, (1) 

where a is similar to Dirac matrix in spin space, vf is the Fermi velocity, a = ( o x .o y , a z ) is the Pauli matrix and the 
Fermi energy is given by Ep. The third term indicates the exchange Hamiltonian ( H ex = a.h), which is present due 
to the interaction between the local magnetization of the ferromagnet and the surface Dirac fermions, where h is the 
exchange energy vector. The Hamiltonian in (1) is not in the form of a usual graphene Hamiltonian [13]. But here 
we have considered pristine spin chiral configuration of ferromagnetic graphene. There are other approaches [14] of 
treating the ferromagnetic character in the Hamiltonian as well. But as we are interested in the spin chiral nature we 
can proceed with the Hamiltonian (1). The term H ex is similar to the Zeeman term which appears as a consequence 
of external magnetic field. In mean field theory one can consider this exchange field as an effective magnetic field. 
The last term on the right hand side of the eqn. (1) is the SOC term, with Xq as the SOC strength. V(r) represents 
the total potential in the system which includes external potential due to external electric field ( 14 ), potential due to 
static disorder (V,i) and crystal potential (Vo)- In this context, the charge and spin Hall effect [15-18] in spin chiral 
ferromagnetic graphene is discussed in [5]. The exchange interaction thus induced can be non uniform [6, 7] as well. 
The exchange vector may be a function of coordinate, momentum or time. The authors in [6], have experimentally 
demonstrated that the exchange coupling in case of the deposition of ferromagnetic insulators may be momentum 
dependent. On the other hand the spin lens is described in [7] with the help of space dependent exchange field. The 
non-uniform exchange coupling can produce interesting effects and can provide a window to investigate enormous 
important issues in the area of spin transport. 


III. CONDUCTIVITY IN FERROMAGNETIC GRAPHENE 

The Hamiltonian of the system with momentum dependent exchange field can be written as 

H = VFOt.k + Ep + d.h(k ) + V(r) + Xg\<J x V r V(r)\.k (2) 

Collecting only the dynamical terms, we can rewrite the Hamiltonian as 

H(k) = Eo{k) + Ep + a.(h + m)(k ) = e(k) + a.M(k ), (3) 

where e(fc) = Eo{k) + Ep, m(k) = AcV r V(r) x k and M(k) = (h + fh)(k). We do our analysis for exchange energy 
d.h(k ) < Fermi energy Ep , such that a is equal to the unit matrix in the spin space and the carriers are electron like 
with up and down spin. Here i = 1,2,3 and e(fc) is the kinetic dispersion energy. Also fh{k) = X cV r F(r) x k and 
M{k) = (h + m)(k). 

In ferromagnet with a constant exchange energy a specific type of charge Hall effect is predicted [5] that is generated 
by spin Hall mechanism in the absence of an external magnetic field. Within the semi-classical theory of spin-orbital 
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dynamics of carriers, a longitudinal electric field produces a pure charge transverse current with no polarization of 
spin. 

Using the Kubo method [19], it is possible to derive the Hall conductivity from Hamiltonian (3). Here we consider 
the exchange vector as momentum dependent and our goal is to show that the conductivity is intimately connected 
to the k space Berry curvature, which is quite similar as the semiconducting system. 

Using the Matsubara Greens function technique, it is well known that it is possible to write the conductivity in the 
following form 


Oxy 


'l 

lim —Q X y{uj + id), 

uj — X) UJ 


( 4 ) 


with 


QxyiiVm) = ^ tr{J x (k)G[k, i{u n + v m )\J y (k)G(k, iu n )}- 


( 5 ) 


k,n 


G(k,iui n ) is the single particle Greens function, SI is the system area, Ji(k) s are the current operator. u> n and v m are 
the fermionic and bosonic Matsubara frequencies and can be expressed as Lo n = (2 n + l)//3 and v m = 2to7t// 3. The 
single particle Greens function is given by 


G(k, iu> n ) = (iuj n — H(k)j 


-l 


1 
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j —h 

iu> n — E+(k) iu} n — E_(k) 


( 6 ) 


where M = y Mi(k)M t (k), R± = |[1 + Mi(k)a l /M\ and E±(k) = e(k) ± M. Using the Hamiltonian the current 
operator Ji(k) can be derived as 


uk) = d ^ = d 4^ + d ^m^ 


dk. 


dki 


dk. 


where i,j = x,y are the space indices. Q xy {wm ) can be defined as 


Qxy(iVm) = -jj-j ^2 ^ 2 tr Jx(k)G{k,i(u} n +i/m))Jy(k)G(k,iu} n ) 


s,£=± 


( 7 ) 


( 8 ) 


k,n 


where s, t are band indices and n is the Matsubara index. We can now recast Q xy (iVm) in terms of the Fermi 
distribution function if we employ the Matsubara frequency sum over the indices n , and we can write Q xy {i v m) as 
[ 20 ] 


Q„(n,) = i E E („, - „.)(*), 


s,t=± % 


iv m - E s (k) + E t (k) 


(9) 


where {n t - n s )(k ) = ( n F {E t (k )) - n F (E s (k)), and n F (e t {k )) = e $E t+1 - Here n F (e s (k )) = e g (E3 _^ m)+1 denotes the 
Fermi distribution function for energy E t and E s respectively. The expression of (Jij can now be written as 


i Tij = lim — Qij(ui + id). 

uj — ►O UJ 


(10) 


Inserting (9), in (10) we can write 
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( 11 ) 


where D(k ) = tr[Ji(k)R_(k)Jj(k)R + (k) — Ji(k)R + (k)Jj(k)R-(k)]. Explicitly in terms of the exchange vector we 


can write 




1 dM a (k) OMpjk) n/r , 
2 M dk x dk y 7 ^ 

2ie a/ 3 7 dM a dMp 
M dk, dk„ 75 


a a^a 1 — a a ' 


( 12 ) 


where a,/3, 7 are the space indices. In this simplification, we have used the relation tr(a a a^a 1 — a a a 7 a^) = 4ie a| g 7 . 
We should remind here that M(k) = (h + m)(k ), which depends on the exchange vector as well as the SOC term. 
Thus the Hall conductivity in the x-y plane can be written as 


_ 1 dM a dMp %/f f 

Txv ~ 2M 3 C V £a/?7 dk x dk v Ml ^ + ){k) - 


(13) 


The dependence of the Hall conductance on the exchange held h(k) is evident in (13). Over the first Brillouin zone 
we can write the Hall conductivity in x-y plane as 


1 


= x ^ 


dM dM 


2 C 


> FBZ 


47T 2 


dk x dky 


where c x c y = C. From (14), we can definitely write the conductivity as 


’xy 


87T 2 


'FBZ 


d 2 kfl z (k), 


(14) 


(15) 


where H z (/c) is the Berry curvature in momentum space. One may note here that for V(r) = Vo, where Vo is the 
crystal potential (which is considered to be constant), the conductivity depends only on the inhomogeneity of the 
exchange vector in momentum space. The Hall conductivity is then given by 


'xy 


Fc 


CxCy -d 2 k 


'FBZ 


4-7T 2 


h. 


dh dh 
dk x dky 


(16) 


In this context it is of interest that even in the absence of an electric Held, momentum dependence of the exchange 
vector can induce Hall conductance in ferromagnetic graphene. In [ 6 ], the momentum dependence of exchange coupling 
is experimentally demonstrated. Our next goal is to choose a specific choice for the momentum dependence, which 
may produce some impressive results. 

The Berry curvature for spin Hall systems carries a spin dependence and is equal but opposite for spin up and down 
electrons i.e fij(fe) = —f2l(fc). Thus eqn (15) indicates that we have the spin Hall conductivity as 


& s h = cr 
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whereas the charge conductivity is given by 


a 


C 

xy 


= o- T 

xy 


+ a 


i 

xyi 


(18) 


which is effectively zero. Thus for a ferromagnetic graphene system the non inhomogeneity in the exchange coupling 
can produce a pure transverse spin current with equal number of spin up and spin down electrons in our system. 


IV. SPIN GAUGE AND THE CURVATURE 

The choice of a specific form of the dependence of the exchange coupling h(k ) on the momentum may help us to 
elucidate the effective Berry curvature in the system from a field theoretical analysis. A simple form of the momentum 
dependent coupling as [21] h(k) = h h{k x , —k v ,0), when incorporated in the Hamiltonian (1) gives, 

H = v F a.k + Ep + h(a x k x — <J y k y ) + V(r) + Ag[ct x V r V(r)].fc. (19) 

The second term reminds us about the well known Dresselhaus coupling in semiconductor. Although the mathematical 
structure of this coupling is similar to that of the Dresselhaus coupling, but the physical origin is far different from 
the original one. This new kind of SOC coupling originates from the momentum dependent exchange coupling of the 
ferromagnetic graphene, that has been experimentally demonstrated in [ 6 ]. The SOC coupling appears here is not 
because of the inversion asymmetry. Let us now consider that the potential V(r) is only due to an external electric 
field in the z direction, so that the relevant part of the above Hamiltonian can be written as 

H = v F a.k + h{a x k x - cr y k y ) + aa{<J y k x - a x k y ), (20) 

where we take Vp = K = 1. A straight forward simplification gives, 

H = a.p + ( hk x - aoky) o x + ( a G k x - hk y ) a y = a.(p- A G ), ( 21 ) 

where Af is the spin dependent gauge field induced by the spin orbit coupling and exchange coupling terms and 
ot G = A g E z , with E z = -V z V(r) 

Explicitly the gauge fields acting in the x and y direction are given by 


A x = -i {-hky + a G k x ) a x (22) 

A!y = i {hk x - a G ky) a y (23) 

The SU(2) gauge field Af is spin dependent and non-Abelian by nature. From the field theoretical consideration this 
gauge field can generate a physical field through the curvature term. The physical field generated due to the presence 
of the gauge A G is given by (H = 1) 

n G = n G y = d x A G - d y A G - i [A G , A g ) (24) 

It is quite obvious that the z component of the k space curvature only exists and is given by 

Hf(fc) = 2 {{-hk y + a G k x ) {hk x - a G k y ) a z ) (25) 

= ±2 (achk 2 - k x k y {h 2 + a 2 G )) , (26) 

where k 2 = k x + k 2 and ± indicates the two spin species of the carriers. The expression clearly shows that the 
curvature is different for the up and down electrons. The variation of this curvature with k x and k y is plotted in figure 
(2) (a). Using equations (15) and (26) we can calculate the conductivity in case of a momentum dependent exchange 
coupling as 


7r a G h 

<Txy = ^r' 


(27) 


where a is the lattice spacing. Thus the conductivity of the ferromagnetic graphene depends not only on the SOC 
coupling parameter but also on the exchange parameter as well. 
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FIG. 1: (Color online) Variation of curvature with momentum. 

V. REAL SPACE BERRY CURVATURE FOR IN-HOMOGENEOUS EXCHANGE COUPLING 

In the previous section we have considered a specific configuration of the exchange vector when it is momentum 
dependent. Now our goal is to observe the consequences when the exchange field is coordinate dependent, which is 
already employed in the study of spin lens configuration [7, 8]. For an non-uniform exchange vector the Hamiltonian 
for ferromagnetic graphene (without any external parameters)can be written as 

H = Vpd.k + a.h(r) (28) 

For the sake of simplicity the spin orbit coupling (SOC) term is also not considered. Actually, our motivation is to 
investigate the role of non-uniform exchange vector on real space berry curvature. Variation of the exchange field can 
induce spin Berry gauge as follows 

= -i (tU(r)|^-| Urtr)} = (29) 

which can be rewritten as 

AiHh) = (tl,h\-^\U,h) (30) 

where a = i,j, k are the space indices. A^(h) is the exchange field dependent Berry gauge field appearing due to the 
inhomogeneity of the exchange field vector. 

A physical field can be generated due to the presence of the gauge A^(h). We only have the ^component of Berry 
curvature as 


n c {r) = 


dAl 1 dJ& _dh a dh b dA? dAA 


dx dy dx dy dh a dhb 
For further analysis, in a standard notation of unit vector we write the exchange vector as 

h = h(sindcos4 >, sin6sin<f >, cosd) 


(31) 


(32) 


where 9 is the polar angle and <\> is the azimuthal angle. 

In the adiabatic approximation, the carriers remain in the same spin eignestates for h(r \)and h{r 2 ), and the flipping 
between states is forbidden. Thus we can write the adiabatic gauge as 


A ld(h) = ±^-(1-cos9)W h (/) 


1 


(33) 


where ± denotes t eigenstate parallel (l anti- parallel) to h(f). Thus the r space Berry curvature is given by, 

^ /-a dha dhjj , h c . 

n ‘ (i1 = -& aT "* 1 S 5 *' 


( 34 ) 
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In the adiabatic approximation, for the coordinate dependent exchange vector a physical field is generated which is 
the well known Berry curvature. This expression of Berry curvature shows that it is the source of monopole, if they 
exists. 

Now for a simplest dependence of the exchange coupling i.e h = ( x , y 1 z) we can write the final form of Berry 
curvature as 


Q = ± 


1 h 

2p’ 


(35) 


which is of nothing but the field generated due to magnetic monopole at the origin. But in general eqn. (31) can be 
written as 


0 , 1 r (dh dh 


(36) 


This is the expression of Berry curvature when we have the exchange field as space dependent. Following [22] we 
can write the Drude Hall conductivity as 


'xy 




(37) 


where n is the electron density and r is the mean free time. Thus with the help of non-uniformity of the exchange 
field we can successfully give a unified picture of conductivity for a SOC Hamiltonian. 


VI. ANTI-DAMPING SPIN ORBIT TORQUE 

Apart from the above results we can argue that besides momentum and space dependence, the exchange field 
can be time dependent as well. In that case we can describe the antidamping spin orbit torque [23] in spin chiral 
configuration of ferromagnetic graphene. This can be an interesting issue regarding construction of magnetic memory 
based devices. To calculate the spin orbit torque we can start with equation (1) but with a time dependent exchange 
field. The torque should carry other contributions from divergence of spin current [24]. But here we are restricting 
ourself to the anti-damping torque only due to the exchange coupling term. We can split the Hamiltonian into three 
parts as 


H = Ho + H ex + H u 


(38) 


where Hq = V(r) = eE.r, and H ex = Ja.h(t), and Hi = njrd.fc — ga.Bsoi(k), where Bsoiik) is the magnetic field 
appearing in the rest frame of the electron due to the spin orbit coupling and J is the exchange coupling parameter. 
Let us consider the electric field in the x direction and the exchange held acts in the opposite to the external electric 
held i.e —h || E. The exchange held induces an additional magnetic held B ex = (2 Jh, 0, 0) in the equilibrium. If we 
consider B ex Bsoi, then we can write the component of spin in the out of plane z direction as [ 11 ] 


s z 


a G eE x 
T 2J 2 h 2 


(39) 


From this equation it is clear that the out of plane spin vector depends not only on the SOC parameter but also on 
the exchange held as well. Also it depends on the direction of the exchange held with respect to the external electric 
held. It is maximum when h and E are anti parallel and zero for parallel orientation. In this generalized picture 
equation (39) can be written as 


Sz,h 


a G eE x 
T 2J 2 h 2 


COsOh'E- 


Integrating s Zt h over all occupied states, the total non-equilibrium spin polarization can be obtained as 


(40) 


S z — 2 gJhs z ^hi 


(41) 


where g is the density of states. This non-equilibrium out of plane spin polarization produces a held in the out of 
plane direction which exerts a torque on the in plane exchange held. The expression of the torque can be written as 


dh 

dt 


— (hx S z z) ~ h x 
h 


E x z) x h 


(42) 





The spin orbit torque can thus be obtained in terms of the exchange field. This can be studied in future. This aspect 
of magnetized graphene may play important role for using magnetized graphene as a future spintronic material. 
Furthermore, following [11] one can argue that the origin of this torque is actually the Berry curvature. The final 
expression of the out of plane spin polarization (non-equilibrium), which is basically responsible for the antidamping 
spin orbit torque can be written as 


S z = ^^2(n+- n~)n z , (43) 

k 

where n± is the Fermi distribution function, V is the volume and Sl z is the z component of the Berry curvature. 
It is hoped that in future this preliminary analysis on the spin orbit torque will be extensively studied and provide 
some impetus to work on the future spintronic memory based devices with ferromagnetic graphene. The beauty of 
this analysis is that it helps us to feel the importance of Berry curvature in analyzing the anti-damping spin orbit 
torque(SOT) in ferromagnetic graphene. One can simply argue that the origin of this SOT is nothing but the Berry 
curvature. 


VII. DISCUSSION 

Induced magnetism in graphene is a gateway of enormous new physics and applications. By depositing EuO on 
top of the graphene one can control the exchange coupling through the molecular beam epitaxy (MBE) method. In 
this paper taking into account the non-uniform exchange field we have theoretically derived the expression of Berry 
curvature in the gauge theoretical approach. For a comprehensive review on gauge theory in spintronics see [25]. 
The non-uniformity might be due to coordinate dependence or momentum dependence. Both cases are studied here. 
With the help of Kubo formula and the Matsubara Greens function techniques we have analytically calculated the 
expression of Hall conductivity in case of a spin chiral ferromagnetic graphene. In the final form of the conductivity 
we clearly can visualize the dependence of the exchange field on Hall conductivity. Consideration of a particular form 
for the non-uniformity we can achieve a synthetic SOC field in the system, which has recent potential applications. 
Furthermore, the time dependence of the exchange vector helps us to explore a new arena of the physics behind 
the spin orbit torque in ferromagnetic graphene. The underlying origin of this torque is nothing but the the Berry 
curvature. In this context, it is of worth mentioning that there would be other components of the torque [24] due to 
the divergence of spin current. But here we are focusing on the torque due to the exchange coupling only. 

The experimental demonstration [6] of the dependence of the exchange vector on the momentum inspired us to 
make the present analysis. Our gauge theoretical analysis shows the importance of the non uniform exchange vector 
in the study of various aspects of spin related issues in spin chiral ferromagnetic graphene . 
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